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Meson Correlation Functions in the e- Regime * 

T. Chiarappa , W. Bietenholz ^, K. Jansen K.-I. Nagai and S. Shcheredin 

NIC/DESY Zeuthcn, Platanenallee 6, D-15738 Zeuthen, Germany 
^ Institut fiir Physik, Humboldt Universitat zu Berlin, Newtonstr. 15, D-12489 Berlin, Germany 

We present a numerical pilot study of the meson correlation functions in the e-regime of chiral perturbation 
theory (xPT). Based on simulations with overlap fermions we measured the axial and pseudo-scalar correlation 
functions, and we discuss the implications for the leading low energy constants in the chiral Lagrangian. 



1. THEORETICAL BACKGROUND 

Ginsparg- Wilson fermions [1,2] obey an exact, 
lattice modified chiral symmetry [3]. Therefore 
they have exact zero modes, which provides to- 
gether with the Index Theorem a conceptually 
clean definition of the topological charge [2] . 

In particular the overlap operator [4] repre- 
sents a relatively simple solution to the Ginsparg- 
Wilson relation. It allows us to tackle now a num- 
ber of physically interesting issues, which were 
inaccessible to numerical studies before. In par- 
ticular, there is hope for simulations to penetrate 
the regime of light mesons, which is described by 
xPT. The latter involves low energy constants as 
free parameters, which can in principle be deter- 
mined by QCD simulations in that regime. To 
the lowest order, the effective chiral Lagrangian 
of xPT takes the form 

Cosp] = ^Tvid^Ud^U] 

- ^Tr 
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(1) 



where U{x) S SU{Nf), Nf is the number of fla- 
vors, M is the (diagonal) quark mass matrix and 
9 the vacuum angle. The pion decay constant 
and the scalar condensate S appear here as the 
leading low energy constants. 

From the practical point of view, the e-regime 
of xPT looks particularly attractive because it 
deals with small volumes [5] . This regime is char- 
acterized by the hierarchy relation 
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so that the pion correlation length is much 
larger than the linear size of the system L. More- 
over, Aqcd can be related to the cutoff of the ef- 
fective theory, Aqcd ^ ^t^Ft^- One identifies the 
matrix field asU = exp{i\/2£_ / Fy^) , where the field 
^(x) describes the light mesons. The e-expansion 
then relies on a non-perturbative treatment of the 
zero-mode, whereas the excitations (the non-zero 
modes of £,{x)) do fit into the volume and can 
be evaluated perturbatively. In the e-regime, the 
topological charge ly of the gauge field plays an 
important role [6]; observables should be mea- 
sured at fixed values of Note that the values 
of and S in the e-regime are the same as in 
the physical situation of an infinite volume. 

Hence the main virtues of Ginsparg- Wilson 
fermions — no additive mass renormalization and 
a clean definition of the topological charge — ren- 
der this lattice fermion formulation ideally suited 
for simulations in the e-regime. However, they 
are numerically so demanding that for the time 
being only quenched QCD simulations are feasi- 
ble. We performed such simulations with overlap 
fermions and the Wilson gauge action at /3 = 6 
on lattices of size 10'^ x 24 and 12"^. The mass of 
the Wilson kernel in the overlap operator was set 
to —1.4, which is optimal for locality [7]. 

Let us mention that Random Matrix Theory 
applied to QCD yields predictions for the low ly- 
ing eigenvalues (EVs) iX of the Dirac operator [8]. 
Simulations in the e-regime confirm these predic- 
tions, at least for the first non-zero EV, if the 
physical volume is not too small, L > 1.2 fm [9].^ 



lower bound for L has also been obtained from the 
condition for the e parameter to be sufficiently small [10] . 
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Fig. 1 shows the probabihty distribution p^^^ of 
this first EV in various topological sectors, as a 
function of the dimensionless variable z = XYIV . 
In particular we see that in the neutral sector 
there is a significant probability pf^ for very small 
EVs. In turned out that this property makes 
measurements very hard, i.e. they would require 
a huge statistics [11]. However, the situation is 
clearly better at non-trivial topology. Therefore 
we concentrate here on the sector \v\ = 1. 
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Figure 1. The probability distribution for the first 
non-zero EV of the massless Dirac operator, ac- 
cording to Random Matrix Theory. 

We measured various types of meson correla- 
tion functions at momentum p — Q and Euclidean 
time t. The corresponding formulae predicted by 
quenched xPT in a volume x T are given in 
Refs. [12]. They were obtained consistently from 
two methods denoted as "rephca" [13] and "su- 
persymmetric" [14]. The vector correlation func- 
tion vanishes identically. The simplest non-trivial 
case is the bare axial correlation function, 

(A^(t)A^(O)), = :^ + 2mS|,|(z)T./ii(T) (3) 



2/il(T) 



' - r + 1/6 , T^t/T , z = mT.V 
I^{z)K^{z) + I^+i{z)K^^i{z) \ + ■ 



where Ii, and Ki, are modified Bessel functions 
and m is the bare quark mass. 

To the first order the pseudo-scalar correlation 
function reads 
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h2{r) = -—[r^-2i 
Si-'°°P(z) = E.(z)- 
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mSeff(P^)F, (4) 



In this case, and also in the scalar correlation 
function, new parameters of the effective chiral 
Lagrangian have to be taken into account, which 
are specific for the quenched approximation: toq 
and a are a scalar mass and a kinetic coupling, 
respectively. For details and for the form of Seff i 
c*^'' and b^_:^\z) we refer to Refs. [12]. 

2. NUMERICAL RESULTS 

We first consider the axial correlation function 
because it only involves the constants E and i^^. 
On the 10^ x 24 lattice the spatial extent is below 

1 fm, and our data show clearly that in such a 
small volume xPT is not applicable, as expected. 

Therefore we now present results from the 
12** lattice, which corresponds to a volume of 
(1.12 fm)'*, i.e. close to the point where the ap- 
plicability of xPT should set in. Our quark mass 
amounts to m = 21.3 MeV, and our statistics 
involves 78 configurations with v = ±.1. 

Fig. 2 (on the left) shows our data for the axial 
correlator. The curve corresponding to eq. (3) can 
be fitted well over some interval that excludes the 
points near the boundary. This fixes -Fn-, which 
enters in an additive constant. On the right of 
Fig. 2 we plot the resulting F^r as a function of 
the number tf oi t values (around T/2) that we 
include in the fit. We find a decent plateau, which 
suggests — up to renormalization — a value of 

= (86.7 ±4.0) MeV . (5) 

On the other hand, S cannot be easily extracted 
from these fits. It is related to the curvature in 
the minimum, but if we vary S over a wide range 
the fit function is hardly affected. In fact. Fig. 

2 (left) shows fits to i/ = 7 points for S vary- 
ing from to (250 MeV)'^, but those curves can 
hardly be distinguished. 

As an attempt to proceed also in that respect, 
we show in Fig. 3 (above) the results for the ax- 
ial and pseudo-scalar correlation functions, which 
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Figure 2. The axial correlation function (left) and 
the values of F^y obtained from fits overtf points 
around the center t — 6 (right). 

are fitted simultaneously. Here the fits involve the 
parameters a, mo and (introduced in eq. 

(4)), in addition to and S. These five param- 
eter fits lead to decent plateaux simultaneously 
for and 1]}^'°°'', as Fig. 3 (below) illustrates, 
j^i-ioop (jgpgjjjjg |-]^g other four parameters, 
and it represents a one loop approximation to S. 
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Figure 3. Above: measured data and fits for 
the axial and pseudo-scalar correlation functions. 
Below: The values Fj^ and 5]1~'°°p obtained from 
five parameter fits. Here tj is the sum of points 
considered in the simultaneous fits for both corre- 
lation functions. 



3. CONCLUSIONS 

We presented results of a pilot study of me- 
son correlation functions in the e-regime, based on 
quenched QCD with overlap fermions. To avoid 
trouble with very low EVs we worked in the sec- 
tor — 1. Although the statistics is still rather 
modest, we recognize some qualitative features. 
If the physical volume is sufficiently large, the 
data can be fitted to the functions predicted by 
quenched xPT. These fits allow for a good deter- 
mination of Ft^ , whereas the evaluation of E is far 
more difficult. 
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